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ABSTRACT 

Tidal friction is thought to be important in determining the long-term spin-orbit 
evolution of short-period extrasolar p lanetary systems. Us ing a simple model of the 
orbit-averaged effects of tidal friction ( |Eggleton et al.|1998 l, we study the evolution of 



close-in planets on inclined orbits, due to tides. We analyse the effects of the inclusion 
of stellar magnetic braking by performing a phase-plane analysis of a simplified system 
of equations, including the braking torque. The inclusion of magnetic braking is found 
to be important, and its neglect can result in a very different system history. We then 
present the results of numerical integrations of the tidal evolution equations, where 
we find that it is essential to consider coupled evolution of the orbital and rotational 
elements, including dissipation in both the star and planet, to accurately model the 
evolution. 

The main result of our integrations is that for typical Hot Jupiters, tidal friction 
aligns the stellar spin with the orbit on a similar time as it causes the orbit to de- 
cay. This tells us that if a planet is observed to be aligned, then it probably formed 
coplanar. This reinforces the importance of Rossiter-McLaughlin effect observations 
in determining the degree of spin-orbit alignment in transiting systems. 

We apply these results to the only observed system with a spin-orbit misalign- 
ment, XO-3, and constrain the efficiency of tidal dissipation (i.e. the modified tidal 
quality factors Q 1 ) in both the star and the planet in this system. Using a model in 
which incrtial waves are excited by tidal forcing in the outer convective envelope and 
dissipated by turbulent viscosity, we calculate Q' for a range of F-star models, and 
find it to vary considerably within this class of stars. This means that using a single 
Q' ', and assuming that it applies to all stars, is probably incorrect. In addition, we 
propose an explanation for the survival of two of the planets on the tightest orbits, 
WASP-12 b and OGLE-TR-56 b, in terms of weak dissipation in the star, as a result 
of their internal structures and slow rotation periods. 

Key words: planetary systems - stars: rotation - celestial mechanics - binaries: close 
- stars: XO-3 - planetary systems: XO-3 



1 INTRODUCTION 

Since the discovery of the first extrasolar planet around a 



solar-type star (Mayor & Queloz|1995 1, observers have now 



jzfij 
£1< 



detected more than 300 planetdj around stars outside the 
solar system. Many of these planets have roughly Jovian 
masses and orbit their host stars in orbits with semi-major 
axes less than 0.1 AU, the so-called "Hot Jupiters" (HJs). In 



both of the giant planet formation scenarios (see Papaloizou 
|& Terquem|2006|for a comprehensive review), core accretion 



and gravitational instability, it is difficult to produce HJs in 
situ. Close-in planets are likely to form in colder regions of 
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1 see http://exoplanet.eu/ for the latest updates 



the protoplanetary disc, much further out (a ~ several AU), 
before a migratory process brings the planet in towards the 
star and to its present location ( Lin et al.|1996 1. 

The formation of systems of giant planets can be 
thought of as occurring in two oversimplified stages ( |Juric fc| 



Tremaine 2008 1 . During stage 1 the cores of the giant planets 



are formed, they accrete gas and undergo migration, driven 
by the dynamical interaction between the planets and the 
gaseous protoplanetary disc (see Papaloizou et al.|2007 for a 
recent review). This stage lasts a few Myr until the gas dis- 
sipates, by which time a population of gas giants may exist. 
If these form sufficiently closely packed then stage 2 follows 
(some evidence in favour of such packing in multiple-planet 
systems is given by [Barnes &: Gree nbcrg 2007[|. This stage 
lasts from when the disc has dissipated until the present, 
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and primarily involves gravitational interactions and colli- 



sions between the planets. Recent studies of stage 2 (Juric 



fc Tremaine|2008||Chatterjee et al.| 2008 For d fc Rasio|2008 1 
have shown that this is a chaotic era, in which planet-planet 
scatterings force the ejection of all but a few (~ 2 — 3) plan- 
ets from the system, in a period of large-scale dynamical 
instability lasting < 10 8 yr. This mechanism can excite the 
eccentricities of the planets to levels required to explain ob- 
servations. 

Planet-planet scatterings tend also to excite the incli- 
nations of the planets with respect to the initial symme- 



try plane of the system (Juric & Tremainc 2008), though 



this has been found to be less efficient than the excitation 
of eccentricity. This potentially leads to observable conse- 
quences via the Rossiter-McLaughlin (RM) effect ( Rossiter 
1924} |McLaughlin|1924| |Winn et alp005 I. 

The RM effect is a spectral distortion of the radial ve- 
locity data that results from the planet occulting a spot 
on the rotating surface of the star as it transits the stel- 
lar disc. High-precision radial velocity data during a transit 
allow a determination of the angle (A) between the sky- 
projected angular momentum vectors of the planetary or- 
bit (h) and stellar spin (fi), through the RM effect. This 
measured value A is not necessarily the same as the inclina- 
tion (or stellar obliquity) i, which is the angle between the 
equatorial plane of the star and orbital plane of the planet 
(defined by cosi = f2 • h), since A is just a sky-projection 
of this angle. Nevertheless, A gives a lower bound on the 
angle between these two vectors, and they are related by 
cos i — cos h cos I p + sin 7* sin I p cos A, where I p and 7* are 
the angles of inclination of the planetary orbital plane and 
stellar equatorial plane, to the plane of the sky. For a tran- 
sit, the orbit must be close to edge-on, so I p ~ 90°, giving 
cos i ~ sin 7* cos A. 

The RM effect has now been used to measure the degree 
of spin-orbit alignment in 11 systems ( |Winn|[2008[ ), though 
this number is expected to grow rapidly over the next few 
years. These systems are currently all consistent with A be- 
ing zero, with the exception of XO-3 b ( jHebrard et al.|2008[ ), 
which is discussed in §8. 1| below. 

It is possible for HJs orbiting a host star which has a 
distant and inclined stellar companion, or massive inclined 
outer planetary companion, to undergo another type of mi- 
gration. This is Kozai migration ( |Wu fc Murray|2003[ ). The 
presence of such an outer companion to an exoplanet host 
star could cause Kozai oscillations, which produce periods 
of extreme eccentricity in the planetary orbit, if various 



conditions arc satisfied (e.g. see section 1.2 of Fabrycky fc 
Tremaine 20071. The subsequent tidal dissipation that oc- 



curs during the periods of small pericentre distance leads to 
gradual inward migration of the planet. It has even been pro- 
posed that a combination of planet-planet scattering, tidal 
circularisation and the Kozai mechanism using outer plan- 



ets, can produce HJs around single stars (Nagasawa et al 



2008). HJs produced from these processes generally have 



their orbital angular momentum vector misaligned with re- 
spect to the stellar spin axis by large angles - occasionally 



larger than 90° ( Fabrycky fc Tremaine|2007 Nagasawa et al. 
[20M] ). 

Misaligned orbits are not predicted from stage 1 alone, 
so if A is measured to be appreciably nonzero in enough sys- 
tems, then it could be seen as evidence for planet-planet 



scattering or Kozai migration. This is because gas-disc mi- 
gration does not seem able to excite orbital inclination 
( Lubow fc Ogilvie|2001||Cresswell et al.|2007| >. Alternatively, 
if observed planets are all found with A consistent with zero, 
this could rule out planet-planet scattering or Kozai migra- 
tion as being of any importance. 

One important consideration is that at such close prox- 
imity to their parent stars, strong tidal interactions between 
the star and planet are expected to cause significant long- 
term spin-orbit evolution, including changes to the value of 
A (actually the true spin-orbit misalignment angle i) over 
time. If tides can change A since the time of formation, then 
we may have difficulty in distinguishing migration caused by 
planet-planet scattering and Kozai oscillations, from gas- 
disc migration. The tidal evolution of such inclined orbits 
must therefore be considered an important goal in planetary 
evolution studies. In this paper we approach the problem of 
studying the effects of tidal friction on such inclined orbits. 



2 PREVIOUS WORK ON TIDAL EVOLUTION 
RELEVANT TO CLOSE-IN PLANETS 

There has been much recent interest into the subject of tidal 
evolution, and our work can be considered to follow on from 
some of the more recent developments in the field. |Witte fc] 
|Savonije| ([2002 ) studied the effects of stellar spin-down on 
tidal evolution. They developed a sophisticated model in- 
cluding the effects of stellar evolution, resonance locking and 
magnetic braking. They studied dynamical tides in solar- 
type stars, and found that the spin-down due to magnetic 
braking causes resonance locking to become more intense, 
which leads to faster orbital decay of a planet. 



Dobbs-Dixon et al. (2004 1 also studied the effects of 



stellar spin-down, though they considered equilibrium tides, 
and we adopt the same model of tidal friction. They stud- 
ied tidal eccentricity evolution, and proposed an explanation 
for the coexistence of both circular and eccentric orbits for 
the planets in the period range 7-21 days, the so-called bor- 
derline planets, as the result of the variation in spin-down 
rates of young stars. Planets with orbital periods less than 
6 days are sufficiently close to their host stars for tidal dis- 
sipation in the planet to be mostly able to account for their 
negligible eccentricities. Those with periods longer than 21 
days are negligibly affected by tides and have been observed 
with a range of eccentricities, so the borderline planets are 
referred to as such since they lie in the crossover between 
these two regimes. Both of these works highlight the impor- 
tance of studying the effects of stellar spin-down on tidal 
evolution. 

Recently |Pont| ( |2008| > discussed empirical evidence for 
tidal spin-up of exoplanet host stars, and found indications 
of such a process occurring in the present sample of transit- 
ing planets. He also proposed that the mass-period relation 
of close-in planets could be accounted for by tidal transfer 
of angular momentum from the orbit to the spin of the star. 
We must note, however, that angular momentum losses from 
the system due to magnetic braking are likely to change this 
picture considerably. In particular, the attainment of a spin- 
orbit synchronous state becomes questionable when angular 
momentum is being removed from the system - this is dis- 
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cussed in ^5]below. We note that his conclusions also depend 
on assumptions regarding the stellar Q' - see Ej3. 1 



|T997||Papaloizou fc Savonije|1997||Ogilvie fc Lin|2004[ here 



Jackson et al.l (|2008a| found that simple timescale con- 



siderations of circularisation may not accurately represent 
the true evolution from considering coupled evolution of the 
eccentricity and semi-major axis. In addition, they found 
that it is inaccurate to neglect the combined effects of the 
stellar and planetary tides in computing the evolution. 

Here we aim to study the accuracy of the simple 
timescale estimates for tidal evolution, when coupled evo- 
lution of the orbital and rotational elements is considered, 
in a more general model of the long-term effects of tidal fric- 



tion than Jackson et al. ( 2008a| consider. We include stellar 



spin-down, and study its effects in a simplified system. In 
particular, we investigate the tidal evolution of inclination, 
since this has not been done in previous studies. 



3 MODEL OF TIDAL FRICTION ADOPTED 

3.1 General introduction to tidal friction 

The tidal interaction between two orbiting bodies acts to 
continually change the orbital and rotational system param- 
eters, and continually dissipates energy. Ultimately - in the 
absence of angular momentum loss from the system - ei- 
ther an equilibrium state is asymptotically approached, or 
the two bodies spiral towards each other at an accelerating 
rate, and eventually collide. The equilibrium state is char- 
acterised by coplanarity (the equatorial planes of the bodies 
coincide with the orbital plane), circularity and corotation 
(rotational frequencies of each body match the orbital fre- 
quency) ( jHut|1980 |. 

The efficiency of tidal dissipation in a body is often 
parametrised by a dimensionless quality factor Q, which re- 
flects the fact the body undergoes a forced oscillation and 
dissipates a small fraction of the associated energy during 
each oscillation period. This is analogous to the quality fac- 



tor iu a forced, damped harmonic oscillator (Murray & Der- 
mott|1999 1, and is defined by 



Q = 2-kEq & -Edt 



where Eq is the maximum energy stored in an oscillation 
and the integral represents the energy dissipated over one 
cycle. This is related to the time lag r in the response of 
the body to tidal forcing of frequency ui by Q" 1 — ujt, when 
Q>1. We find it convenient to define Q' = I?, where k is 
the second-order potential Love number of the body, since 
this combination always appears together in the evolution- 
ary equations. Q' reduces to Q for a homogeneous fluid body, 
where k = §. 

The problem of determining the efficiency of tidal dis- 
sipation, and therefore quantifying the evolution of the sys- 
tem, amounts to calculating Q' factors for each body. Q' is 
a function of the tidal frequency (and possibly the ampli- 
tude of the tidal disturbance), being the result of complex 
dissipative processes in each body ( |Zahn|[2"b08j ). For rotat- 
ing fluid bodies, such as giant planets and stars, calculations 
of the excitation and dissipation of internal waves have in- 
dicated that Q' varies in a complicated way with the tidal 
frequency (see Savonije et al.|[l995 Savonije & Papaloizou 



after OL04; |Ogilvie fc Lin|2007| hereafter OL07). These cal- 
culations rely on a variety of approximations and involve 
uncertainties, particularly regarding the interaction of the 
waves with convection. 

Typically assumed values of Q' ~ 10 6 for stars are 
roughly consistent with observational data regarding the cir- 
cularisation periods of binary stars (OL07). In addition, the 
magnitude of Q for HJs is often assumed to be similar to 
that for Jupiter, w hich has been inferred to be in the range 
6 x 10 4 - 2 x 10 6 (lYoder & Peale||l98l|). This estimate is 



based on a model of the tidal origin of the Laplace reso- 
nance among the Galilean satellites; however, it has been 
argued that even if the origin of the resonance is primor- 



dial, the ave rage Q cannot be far from these bounds ( Peale 
fc Lee||2002[ ) - giving Q' ~ 10 6 (since k ~ 0.38 for Jupiter 



This estimate also appears consistent with the work of |Jack-| 



son et al. (2008a I, who found that one can reproduce the 



outer planet (o > 0.2 AU) eccentricity distribution from in- 
tegrating the tidal evolution equations backwards in time 
for the observed close-in planets (a < 0.2 AU) quite well if 
Q ~ 10 5 5 for stars and Q ~ 10 6 ' 5 for HJs. However, this 
stellar Q' is difficult to reconcile with the existence of the 



planets on the tightest orbits, such as WASP- 12 b (Hebb 
|et al.|[2008] | and OGLE-TR-56 b ( |Sasselov|[2003l ), since it 
would imply that the inspiral time for these planets would 
be much less than the age of the system. That several plan- 
ets have been found with similarly short periods makes this 
seem unlikely on probabilistic grounds. Additionally, |Jack-| 



son et al. (2008a I assume that Q' is the same for all exo- 
planet host stars; this may be an oversimplification, as is 
discussed in §8. 2| below. 

It would seem more plausible that the relevant stellar 
Q' is in fact higher than this, and this is supported by the 
theoretical work of OL07. They also propose that the effi- 
ciency of tidal dissipation in solar-type stars may be dif- 
ferent when the orbiting companion is a ~ Mq star in a 
close binary, than when its orbiting companion is a close-in 
gas giant planet. This could be a result of differences in the 
tidal and spin frequencies in the two situations, which may 
allow the excitation of inertial waves in the convective en- 
velope of the star in the former case but usually not in the 
latter. It could also be a result of the excitation of internal 
inertia-gravity modes at the boundary between the convec- 
tive and radiative zones, and their resulting dissipation in 
the radiative core of solar-type stars. If these waves are of 
small amplitude, then they are unlikely to achieve sufficient 
nonlinearity to prevent coherent reflection from the centre 
of the star, hence they form global modes, and only weakly 
contribute to the dissipation ( Goodman fc Dickson|1998 i. If 
these waves achieve sufficient nonlinearity, then they could 
overturn the local entropy stratification. This may result in 
wave breaking, which could dissipate energy in the wave, or 
in incoherent reflection of the waves from the centre of the 
star, preventing the formation of global modes, and enhanc- 
ing the dissipation. Simple estimates of when these waves 
become nonlinear indicate that the waves excited by HJs 
marginally achieve sufficient nonlinearity to disrupt the re- 
flection of the waves. This is in contrast to the close binary 
circularisation problem, in which these waves are likely al- 
ways to achieve sufficient nonlinearity. This may result in 
much higher Q' values relevant to the survival of the very 
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close-in planets, potentially explaining the discrepancy be- 
tween the circularisation of binary stars and the survival of 
the very close-in planets. 



3.2 Model adopted — see Appendix [A] for details 

In light of the uncertainties involved in calculating Q' , and 
the difficulty of calculating the evolution when Q' is a com- 
plicated function of ui, we adopt a simplified model, based 
on a frequency-independent lag time - though see discus- 
sion in Appendix [A] We adopt the model of Egglcton ct al. 
( 1998), which is based on the equilibrium tide model oflHut 



( 1981 1. In this formulation, we calculate the evolution of the 



specific angular momentum of the planetary orbit 
h = r x r = na 2 vl-e 2 h, 



together with its eccentricity vector e, and the stellar and 
planetary spin vectors fii and CI2. The eccentricity vector 
has the magnitude of the eccentricity, and points in the di- 
rection of periastron, and is defined by 

rxh 

e = 79 r, 

where m\2 = vn\ +rri2 is the sum of the stellar and planetary 
masses. 

Both h and e are conserved for an unperturbed Keple- 
rian orbit; therefore under weak external perturbations their 
components vary slowly compared with the orbital period. 
This allows averaging of the effects of the tidal perturbation 
over a Keplerian orbit, resulting in a set of secular evolution 
equations for the rotational and orbital elements. 

This formulation is beneficial because it can treat arbi- 
trary orbital eccentricities and stellar and planetary obliq- 
uities, unlike other models which are only valid to a given 
order in the eccentricity, or for small (or zero) orbital incli- 
nations ( Goldreich fc Soter|1966 Jackson et al.|[2008a Hut 
19811. Using the secular evolution equations allows us to 



perform integrations quickly that represent dynamical evo- 
lution over billions of years. The full set of equations is pre- 
sented in Appendix [A] in a form which is straightforward to 
numerically integrate. These equations have been written in 
such a way as to eliminate references to the basis vectors 
chosen in their representation, since the eccentricity basis 
vector is undefined for a circular orbit. In this form, the 
equations are regular at e = 0, unlike those in |Eggleton fc| 
Kiseleva-Eggletonl (120011) andlMardling & Lin|d2002l. 



4 MAGNETIC BRAKING 

Observations of solar-type stars have shown that the mean 
stellar rotational velocity decreases with time (Skumanich 
1972), following the relation 51 oc t~ 1//2 , where t is the 
main-sequence age. This is the empirical Skumanich rela- 
tion, and can be interpreted as telling us that solar-type 
stars have been undergoing continuous spin-down since they 
first started on the main sequence. Magnetic braking by a 
magnetised outflowing wind has long been recognised as an 
important mechanism for the removal of angular momen- 
tum from rotating stars (Weber & Davis 19671, and such 
a mechanism seems able to explain most of the observed 



Although the Skumanich law is well established for stars 
with rotational velocities in the range 1-30 kms" 1 , it over- 
estimates the spin-down rates of stars up to t ~ 10 s yrs, 
and thus cannot explain the presence of fast rotators in the 
Pleiades (Ivanova & Taam 2003). As a resolution to this 
problem it has been suggested that the angular momentum 
loss rate for high rotation rates could be reduced, as a result 
of the saturation of the stellar dynamo ( MacGrcgor & Bren 



1991 1, or alternatively due to a reduction in the num- 



ber of open field lines in a complex magnetic field topology 



I Taam & Spruit 19891. These and similar approaches lead 



to modified models of the magnetic braking torque for fast 
rotators, and several such models have been proposed (e.g. 
Ivanova fc Taam||2003| |Holzwarth fc Jardine||2005[ >. Never- 
theless, Barnes (20031 finds that the Skumanich relation is 



remarkably accurate at modelling the spin-down of Sun- 
like stars that are not rapid rotators, so to a first approx- 
imation, a magnetic braking torque based on the empirical 
Skumanich law is best for our purposes. 

Here we include the effects of magnetic braking in the 
tidal evolution equations, through the inclusion of the |Ver-| 
bunt & Zwaan ( 1981 1 braking torque, with the particular 
coefficients of Dobbs-Dixon et al. (20041, as follows (see Ap- 
pendix [a} 



limb — —CXmh ^1 ^1 



(1) 



where a m b = 1-5 x 10~ 14 7 yrs. 7 is a correction factor for 
F-dwarfs, which takes the value 0.1 for an F dwarf, but is 
unity for a G or K dwarf. We can also define a magnetic 



braking timescale r m b 
for the Sun. 



^2 , which is ~ 10 10 yrs 



5 ANALYSIS OF THE EFFECTS OF 
MAGNETIC BRAKING ON TIDAL 
EVOLUTION FOR A SIMPLIFIED SYSTEM 

5.1 Circular, coplanar orbit with magnetic 
braking 

We first study the effects of magnetic braking on a simpli- 
fied system of a circular, coplanar orbit under the influence 
of only the tide that is raised on the star by the planet, 
and magnetic braking. We have neglected the tide in the 
planet here since the moment of inertia of the planet is much 
smaller than that of star and the orbit (i.e. I2 <C Ii ~ m,2a 2 ), 
so to a first aproximation we can neglect the effects of plan- 
etary spin; in any case the planetary spin is expected to 
synchronize rapidly with the orbit. The following set of di- 
mensionless equations can be derived from the full set of 
equations in Appendix [A] 



dt 
dn 
dt 



n 



iff 



= 3n"3" 



(2) 
(3) 



stellar spin- down ( |Barnes|200'3" l 



where we have normalised the stellar spin frequency Qi and 
orbital mean motion n to the orbital frequency at the stel- 
lar surface, together with a factor C 3//4 . C is the ratio of the 
orbital angular momentum of a mass m,2 in an orbit with 
semi-major axis equal to the stellar radius _Ri , to the spin an- 
gular momentum of an equally rapidly rotating star of radius 
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ill, mass mi and dimensionless radius of gyration r g \. The 
reduced mass is u = mi , m2 . C is important for classifying 

the stability of the equilibrium curve Cl = n in the absence 
of magnetic braking, and it can be shown from energy and 
angular momentum considerations that this equilibrium is 

_ 3 

stable if h ^ 3 4 - equivalent to the statement that no more 
than a quarter of the total angular momentum can be in the 



form of spin angular momentum for stability (Hut 1980). We 



have thus defined the following dimensionless quantities: 



Q = fii 



( Rl 



c 



\Gmi2 
( Rl 



C" 



h 

( Gm\2 



mi 



9 
Wi 



m 2 
mi 



CT t, 



A 



= Qmb 



/ Gmi 



\ Rl 



2Qi 
9 



mi 
ma 



C" 



There is only one parameter (A) that completely char- 
acterises the solution in the (n, f2)-plane, and its value may 
be estimated as 



A ~ 100 7 



10 6 



for a Jupiter-mass planet orbiting a Sun-like star undergo- 
ing magnetic braking (with standard a m b and with Q' — 
10 6 ). The size of this term shows that in general magnetic 
braking dominates the stellar spin evolution. Note that in 
the absence of magnetic braking (A = 0), Eqs. [2] and [3] do 
not contain reference to the masses of the star and planet 
or to the tidal Q' of the star. The parameter A, together 
with the initial conditions (no,fio), completely determines 
the evolution. 

We plot some solutions on the (n, fi)-plane in Fig.[l]by 
integrating Eqs. [2] and [3] for various initial conditions. Fig.[l] 
shows two phase portrait plots, which show the general qual- 
itative behaviour of the solutions to Eqs.[2]and|3j for a given 
value of the parameter A. The arrows on each curve show 
the direction of time evolution from the initial state. For 
prograde orbits we restrict ourselves to studying the region, 
^ n ^ 10, ^ ^ 10, in the (n, f2)-plane. This is because 
fii = \jGm\2l Rl corresponds to stellar breakup velocity 
and n ^ \J Gm^/ R\ means that the planet would be orbit- 
ing at, or beneath, the stellar surface. For a HJ with a mass 
of Mj orbiting a star of mass Mq, C ~ 0.01, so h ~ 10 cor- 
responds to an orbital semi-major axis of a ~ 0.01 AU, and 
n ~ 0.1 corresponds to a ~ 0.2 AU, so these plots represent 
the full range of orbits of the HJs. The top plot in Fig. [I] 
is for A = 100, which corresponds to canonical magnetic 
braking for a G/K star (7 = 1) and Q[ of 10 6 . 

In the absence of magnetic braking (A — 0) we recover 
the standard tidal evolution equations for a coplanar, circu- 
lar orbit. These have been well studied in the literature (e.g. 
Counselman 1973 Hut|1981 1. These equations have an equi- 
librium of coplanarity and corotation i.e. i = and Cl — h, 
where the orbital inclination (or stellar obliquity) i is de- 
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Figure 1. (n,f2)-plane with A = 100 for a HJ orbiting a Sun- 
like star. The diagonal dashed line in each plot corresponds to 
corotation (Q = h). Top: magnetic braking spins the star down 
so that the planet finds itself inside corotation, where the sign of 
the tidal torque changes, and planet is subject to tidally induced 
orbital decay. For an initially high n outside corotation tidal fric- 
tion efficiently transfers angular momentum from spin to orbit, 
which pushes the planet outwards. Bottom: Solutions with the 
same initial conditions are plotted with and without magnetic 
braking for a HJ around a solar-type star, with dot-dashed lines 
having A = and solid lines have A = 100. The dot— dashed lines 
are also curves of constant total angular momentum. This shows 
that the inclusion of magnetic braking is extremely important in 
determining the secular evolution of the system, and its absence 
results in a very different evolutionary history unless Q h in 
the initial state. 



fined by cosi = f2i • h and i ^ 0. The system will approach 
this equilibrium if both the spin angular momentum is less 
than a quarter of the total angular momentum, and the total 



angular momentum exceeds some critical value ( Hut||1980 



Greenbcrg 1974J. With no braking, orbits initially outside 
corotation (SI > h) are not subject to tidally induced or- 
bital decay, and asymptotically approach a stable equilib- 
_ 3 

rium Q, — h for h ^ 3 4 . Orbits initially inside corotation 
can evolve in two different ways, depending on the stability 
of the equilibrium state on the solution's closest approach to 
O = n. Iff^>^>0 near corotation, then n < 3~ 4 , and 
the equilibrium state is locally stable (though no such curves 
are plotted in Fig. [T] since they occur only in the far bottom 
left of the plot, near the origin). This is when the corota- 
tion radius moves inwards faster than the orbit shrinks due 
to tidal friction, which can result in a final stable equilib- 
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>a 5- 

4- 




Figure 2. (n, Q)— plane showing the effects of reducing the stel- 
lar mass fraction participating in angular momentum exchange 
with the orbit while the braking rate is unchanged, as has been 
proposed for an F star like r Boo. A = 100, e* = 10~ 2 



So far we have considered the whole star to participate 
in tidal angular momentum exchange with the orbit. For an 
F-dwarf (like r Boo), it has been proposed that only the 
outer convective envelope (of mass fraction ~ e*) partici- 
pates in angular momentum exchange with the orbit (iMarcy 



et al. 1997 Dobbs-Dixon et al. 2004 Donati et al. 2008) 



If the core and envelope of such a star can decouple, then 
tides would only have to spin up the outer layers of the star, 
which would reduce the spin-up time by ~ e*. In this case 
the system could remain in a state with ~ n just inside 
corotation, with the resulting torque on the orbit small. This 
may explain the spin-orbit synchronism of stars such as r 

p004||. Fig. [2] shows 



Boo, as noted by Dobbs-Dixon et al. 



the phase plane for a simplified system in which the moment 
of inertia of the star acted on by tides is reduced by a factor 



10" 



but the braking rate is unchanged i.e. we mul- 



tiply the first term on the right-hand side of Eq. [2] by el 
Note that this may be too simple a model to describe such 
core-envelope decoupling, and we have ignored associated 
changes to the braking rate. 



rium state for the system if the corotation radius "catches 
up" with the planet. On the other hand, orbits inside coro- 
tation for which this condition is not satisfied are subject 
to tidally induced orbital decay, since tidally induced an- 
gular momentum exchange enhances the difference between 
Cl and n, which leads to further orbital evolution, and the 
spiralling in of the planet. This evolution can be seen from 
the dot-dashed lines in the bottom plot in Fig. [I] 

Including magnetic braking (^4 ^ 0) means that Cl — h 
is no longer an equilibrium state, and the total angular mo- 
mentum of the system is not conserved. For an orbit initially 
not subject to spiralling into the star via tidal transfer of an- 
gular momentum from orbit to spin (O ^ h) we see from the 
top of Fig. [T] that magnetic braking will spin the star down 
so that the planet finds itself inside the corotation radius of 
the star. Passing through corotation changes the sign of the 
tidal torque and causes the planet to spiral into the star. 
Note that, if we ignore the age of the system, any bound 
orbit will eventually decay in a finite time since the system 
has no stable equilibrium. The effect of magnetic braking is 
to increase the minimum semi-major axis at which the orbit 
is not subject tidally induced orbital decay over the nuclear 
lifetime of the star. 

This means that an initially rapidly rotating G-type 
star hosting a close-in Jupiter mass companion will lose sig- 
nificant spin angular momentum through magnetic braking 
(over a time ~ r m b). During this stage of spin-down the 
spin frequency of the star may temporarily equal the orbital 
frequency of its close-in planet, but the rate of angular mo- 
mentum loss through magnetic braking will exceed the tidal 
rate of transfer of angular momentum from orbit to spin. 
The stellar spin continues to drop well below synchronism 
until the efficiency of transfer of tidal angular momentum 
from orbit to spin can compensate or overcompensate for 



the braking. If M > 



dh 

di ' 



inside corotation, then tides 



will act to spin up the star, though the timescale for this to 
cause significant spin-up may be much longer than the stel- 
lar lifetime, and this only occurs if the orbit has sufficient 
angular momentum to noticeably spin up the star. Other- 
wise, the planet continues to spiral inwards once it moves 
inside corotation, and Q ~ const. 



5.2 Extending the analysis to inclined circular 
orbits 

We can extend the simplified system of equations analysed 
in the previous section to arbitrary inclination (i) of the 
orbital plane with respect to the equatorial plane of star: 



di 
dh 
di 
di 
di 



cos i (l + cos 2 i) 

2n v J 



A fi 



= 3nTT 



n 



— I cos 1 — n 3 i I 
2n 



(4) 
(5) 
(6) 



For small inclination, Eq. [6] reproduces Eq. 13 from |Hut| 
( 1981 1, with the exception that we have used a constant Q[ 



rather than a constant time-lag in the equations (i.e. replace 



time lag r by 



- , and note that k± is twice the apsidal 



2 1h\ nQ 

motion constant of the star). 

From Eq. [5] the orbit begins to decay if 

17 cos i < h, 



(7) 



which is always satisfied for a retrograde orbit (i ^ 90°). 
This is just a generalisation of the corotation condition Cl = 
n to a non-coplanar orbit. The inclination grows if 



Q. > Q crit — 2h ( cos i — h 3 Q 



(8) 



where we have assumed the quantity in brackets is posi- 
tive i.e. cosi > n 3 f2. This agrees with the condition from 
Hut ( 1981 1 when i ~ 0. For sufficiently close-in orbits that 



tidal friction is important, magnetic braking will rapidly spin 
down the star such that this condition is not satisfied, so we 
can safely conclude that the inclination is not likely to grow 
appreciably by tidal friction. When this condition is not sat- 
isfied, the inclination decays to zero, on a timescale n (see 
next section). Note that i — 180° is an unstable equilibrium 
value of the inclination. 

In the absence of magnetic braking we recover the evo- 



lution considered by Greenberg (19741, so we will now con 
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centrate on the inclusion of magnetic braking. The top figure 
in Fig. [3] plots the (n, f2)-plane for A = 10 for an initial in- 
clination of i = 90° . A smaller value of A is chosen over the 
previous section in order to show the effects of tides more 
clearly, since reducing A is equivalent to reducing Q' . The 
evolution of i is not plotted (and will in general be different 
for each curve), but is found to decay once the stellar spin 
decays sufficiently that Eq. [8] is not satisfied. The bottom 
figure in Fig. [3] shows the effect of increasing the inclination 
in steps to illustrate the behaviour as i is increased, on var- 
ious curves with otherwise the same initial conditions, with 
A — 10. The orbit begins to decay for smaller h, and decays 
at a faster rate as i is increased. This peaks for a perfectly 
retrograde orbit (i = 180°), with anti-parallel spin and or- 
bit, where the rates of change of spin and orbital angular 
frequencies are maximum. Fig. [3] shows that the orbit gen- 
erally begins to decay outside corotation, once Qcosi < n. 
This has implications for the tidal evolution of close-in plan- 
ets on inclined orbits, in that if this condition is satisfied, 
the planet will be undergoing tidally induced orbital decay 
- though the inspiral time may be longer than the expected 
stellar lifetime. 

In this section we have seen that magnetic braking can 
only be reasonably neglected for a coplanar orbit when Q <C 
n. For an inclined orbit, this condition must be generalised 
to S7cosi <C n - obvious from Eq. [5] Neglecting w m b for 
stars for which this condition is not satisfied can result in a 



qualitatively different evolution, as already seen in Fig. 1(b) 
for a coplanar orbit. 



6 TIDAL EVOLUTION TIMESCALES 

It is common practice to interpret the effects of tidal evolu- 
tion in terms of simple timescale estimates. The idea behind 
these is that if the rate of change of a quantity X is ex- 
ponential, then X/X will be a constant, so we can define 
a timescale tx = X/X. If X/X ^ const, then these may 
not accurately represent the evolution. Here we reproduce 
the timescales that can be derived from the equations in 
Appendix [X] 

A tidal inspiral time can be calculated from the equa- 
tion for d, by considering only the effects of the tide raised 
on the star by the planet (not unreasonable since I2 <C Ii ~ 



m-2a ). Here a/a ■ 



-13/2 



7^ const, so a more accurate esti- 



mate of the inspiral time for a circular coplanar orbit is 
2 a 



13 d 



±12.0 Myr 



9i 

10" 



mi 



Mj_ 
rn 2 



1 - 



Here P and P* are the orbital and stellar spin periods, 
respectively. We have already seen in ^5] that it is unreason- 
able to assume that SI is fixed unless Q <C n, due to magnetic 
braking. 

If the orbit is inside corotation, angular momentum will 
be transferred from the orbit to the spin of the star, giving 
a tidal spin-up time of 

fii 13r a 



2 a 



where a = is the ratio of orbital to spin angular mo- 
mentum. For the HJ problem, tq 1 ^ r a since a ~ 0(1), 




<G 5 




Figure 3. Top: (n, Q)— plane for an orbit with an initial i = 90°, 
with A = 10, where the dashed line corresponds to f2 = n. This 
is similar to Fig. |l(a)| for a circular orbit, except that the orbit 
decays once f2 cos i < h, which can occur above the dashed line in 
contrast with the circular case. The evolution of i is not plotted, 
and is in general different for each curve, but is found to decay 
once the stellar spin drops below that given by Eq. [8] Bottom: 
various initial inclinations, with A = 10, where the dashed line 
corresponds to f2 = n. The bottom trajectory has i = 0, and the 
inclination is increased in steps towards the top curve, which has 
i = 180°, to illustrate the behaviour. Note that orbits with larger 
initial i decay for smaller n, once f2cosi < h is satisfied. 



which neglects the spin-down effects of magnetic braking. 
The planetary spin 172 will tend to synchronise much faster, 
since the moment of inertia of the planet is much less than 
that of the orbit (by ~ 10 5 ), and will not be considered 
further i.e. we assume £^2 = n. 

A circularisation time can be derived from the equation 
-lfor e, and is given for a coplanar orbit by 



16.8 Myr 



10V 



mi 



/i(e2) -i8i /2(e2) 



/MA /Ro\ 5 fP_ 
\m 2 ) \ Ri ) Ud 



/l(e 2 )-^/ 2 (e 2 ) 



where we have included both the stellar and planetary tides, 
as these have been shown to both contribute to the tidal 
evolution of e (Ja ckson et al.|2008a l. Note that e/e ~ const 



only if P ~ const and e< 1, where /i,2(e ) ~ 1. The factor 



= 



mi 
m 2 



10- 
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for the HJ problem. 

If the orbital and stellar equatorial planes are mis- 
aligned, then dissipation of the tide raised on the star by 
the planet would align them on a timescale 



dl 



Vioe 



1 - 



P 
2R 



mi \ 




Mq) 











P\ 4 
id 



where we have assumed that the orbit is circular and made 
the small i approximation. We take fio = 5.8 x 10 -6 s -1 , 
which corresponds to a spin period of ~ 12.5 d. 

The validity of these timescales to accurately represent 
the tidal evolution of the orbital and rotational elements 
is an important subject of study, since these timescales are 
commonly applied to observed systems. In a recent paper, 
Jackson et al. (2008a I found that it is essential to consider 



the coupled evolution of e and a in order to accurately model 
the tidal evolution, and that both the stellar and plane- 
tary tides must be considered. They showed that the actual 
change of e over time can be quite different from simple 
circularisation timescale considerations, due to the coupled 
evolution of a. In the following we will consider the validity 
of the spin-orbit alignment timescale to accurately model 
tidal evolution of i. 



7 NUMERICAL INTEGRATIONS OF THE 
FULL EQUATIONS FOR AN INCLINED 
ORBIT 

We perform direct numerical integrations of the equations 
in Appendix |A| with a 4th/5th order Runge-Kutta scheme 
with adaptive stepsize control, using a scheme with Cash- 
Karp coefficients similar to that described in |Press et al.| 
(19921. Our principal aim is to study inclination evolution 



and to determine the accuracy of the spin-orbit alignment 
timescale n for close-in planets. 

We choose a "standard" system of a HJ in orbit around 
an FGK star. We have mi = Mq, mi = Mj, Ri = Rq, R2 = 
Rj, and we modifj^] a, i, e. We choose the initial ratio 
Qi/n = 10, and include magnetic braking in all simula- 
tions, setting 7 = 1 unless stated otherwise. We choose 
Qi — Q'i = 10 6 , and fli/n = 1 unless stated otherwise. The 
dimensionless radii of gyration are chosen to be r^i = 0.076 
and Tg 2 = 0.261, which are values appropriate for polytropic 
stellar and planetary models with respective indices 3 and 
1. 



7.1 Inclusion of magnetic braking and the 

importance of coupled evolution of a and i 

For a prograde orbit (i < 90°) initially outside corotation, 
w m b rapidly spins the star down sufficiently to ensure that 
inclination is not excited through tidal friction (so that Eq.[8] 



2 We take f2i • e = 0, though this need not be assumed. We find 
negligible difference between integrations for which f2i - e = and 
Sii • e + 



is not satisfied), and the inclination begins to decay. Subse- 
quent spin-down moves the orbit-projected corotation ra- 
dius beyond the orbit of the planet (so that ilcosi < n), 
and the resulting tidal inspiral accelerates as the difference 
between Q cos i and n is enhanced. The associated reduction 
in a increases the rate of stellar spin-orbit alignment. Thus 
the inclusion of ti> m b increases the rate of alignment, and 
reduces n from the simple estimate, which ignores magnetic 
braking and coupled a and i evolution. 

The effect of <i> m b on a retrograde orbit [i ^ 90°) is qual- 
itatively different. A retrograde orbit is always subject to 
tidally induced inspiral since Q cosi < n for all i ^ 90°. u m b 
acts to reduce Q, thereby reducing the difference |f2cosi— n\, 
making the tidal torque smaller. This acts to increase the 
timescale for alignment of the stellar spin and orbit, though 
the effect is found to be small. 

The most important effect of including (i> m b is sim- 
ply that of reducing the stellar spin sufficiently so that 
f2cosz < n, where the semi-major axis can then decay 
through tidal friction. As a subsequently decreases, the tidal 
torque increases, resulting in a faster inclination decay rate. 
w m b can only be neglected if |ficosi| <C n. If the orbit is 
already highly inclined, such as i 90°, then the inclusion 
of magnetic braking is not so important, since ilcosi < n 
regardless of the spin rate. 

For an orbit initially at a — 0.05 AU, the simple es- 
timate of the stellar spin-orbit alignment timescale gives 
Ti ~ 2.0 x 10 10 yrs. Fig. [4] shows the evolution of i and 
a for various initial inclinations for an orbit at a = 0.05 AU 
and a = 0.08 AU respectively. The outer orbit a changes 
only slightly over 10 Gyr, and i evolves as expected from 
the simple estimate of n ~ 3.0 x 10 11 yrs. The inner orbit, 
on the other hand, is subject to tidally induced orbital de- 
cay, with an inspiral time of 1 — 3 Gyr, and this reduction 
in a increases the rate of inclination evolution. This high- 
lights the importance of considering coupled evolution of a 
and i, especially for large initial inclinations, where inspiral 
occurs for higher stellar spin rates. The difference between 
the simple estimate of n and the true timescale can be up to 
an order of magnitude different for orbits whose a changes 
appreciably. 

7.2 Inclined and eccentric orbits 



Jackson et al. ( 2008a I already highlighted the importance of 
coupled a and e evolution for a coplanar orbit. We will now 
consider how e might affect i evolution for a non-coplanar 
orbit, which is the subject of this section. 

A nonzero eccentricity reduces the pericentre distance 
r p = a(l — e), which increases the tidal torque over a cir- 
cular orbit, since tidal torque ~ r -6 . Although the planet 
would spend less time near pericentre, the torque there is 
much greater, so dominates the orbit-averaged torque. We 
therefore expect the stellar spin-orbit alignment time to be 
reduced as we increase e. In addition, we expect that an orbit 
at large e would more strongly affect the rate of alignment 
over a circular orbit, than one at large i would over a copla- 
nar orbit, because e reduces r p , whereas i only changes the 
difference (57 cosi — n). The tidal torque ~ r~ D (f2 cos i — n), 
which depends more strongly on r than Q. This behaviour 
can be seen in Figs. |4(d)] and [5(b)| which shows that the ra- 
tios of stellar spin-orbit alignment times for an orbit with 
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Figure 4. Tidal evolution for a circular, inclined orbit at a = 0.05 AU and a = 0.08 AU, with various initial inclinations: i = 
6° , 26° , 45° , 90° , 180° . (a) and (c) show their respective semi-major axis evolutions (with the highest inclination orbit decaying first - the 
bottom curve - and the lowest inclination orbit decaying last - the top curve), and (b) and (d) show the respective inclination evolution 
for these systems. The outer orbit a changes only slightly over 10 Gyr, and i evolves as expected from simple estimates of rj. The inner 
orbit, on the other hand, is tidally shrunk, and this reduction in a means that the true evolution is much faster than the simple timescale 
estimates predict. This highlights the importance of considering coupled evolution of a and i. 



small i and large i is ~ 0(1), whereas the ratios of stellar 
spin-orbit alignment times for an orbit with small e and 
large e can be up to several orders of magnitude. 



Prom Fig. [5] we can compare the simple estimate of 
n ^ 3.0 x 10 11 yrs for an orbit at a = 0.08 AU, with the 
coupled evolution of the orbital and rotational elements from 
integration of the full equations. We see that the simple es- 
timate gives a misleadingly long stellar spin-orbit alignment 
time compared with that obtained from integrating the full 
equations in cases where e is initially non-negligible - as 
a result of the strong functions of the eccentricity in this 
model of tidal friction (see Appendix [Aj . This confirms the 
conclusion of Jackson et al. (2008a I, in that it is essential 
to consider coupled evolution of a and e in order to de- 
termine an accurate system history. We also find that the 
associated changes in semi-major axis strongly affect stellar 
spin-orbit alignment. A marginally better estimate for n 
can be made by replacing the orbital period with the orbital 
period around periastron (equivalent to replacing a by r p ), 
though this is still inadequate since it neglects evolution of 
a. 



7.3 Discussion 

For typical HJs, we find that the stellar spin-orbit align- 
ment time is comparable to the inspiral time i.e. Ti ~ r n . 
This means that if we observe a planet, then its survival 
implies that tides are unlikely to have aligned its orbit. For 
planets on an accelerating inspiral into the star, the rate 
of inclination evolution will have been much lower in the 
past. Therefore if we observe a planet well inside corotation 
(ficos i < n), with a roughly coplanar orbit, we can assume 
that it must have started off similarly coplanar - unless we 
are lucky enough to be observing a planet on its final rapid 
inspiral into the star after it has undergone most of the evo- 
lution, where it is now in a very short-period orbit, close to 
being consumed. 

We expect r< ~ r a when a< 1- which is true for close- 
in terrestrial planets - since $1 can be considered fixed, with 
the inclination changing only due to changes in h. For typical 
values of a ~ 0(1) for HJs, the inclination changes due to 
rotations of both fl and h, so the timescales are not exactly 
the same, but would be expected to be of the same order of 
magnitude. 

|Hut| ( |1981[ ) showed by considering only the tide in the 
star, that the stellar spin-orbit alignment timescale is longer 
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Figure 5. Top left: Semi-major axis evolution for an inclined orbit with initial i = 45° at a = 0.08 AU for various initial e, with 
e = 0, 0.2, 0.4, 0.5, 0.6, 0.7, 0.8. Top right and bottom: Inclination and eccentricity evolution for the same systems. Solutions with the 
smallest initial e have the smallest change in a and i. Solutions with the largest initial e undergo much more rapid tidal evolution (note 
that for e = 0.8, e decays to less than 0.7 within 10 5 yrs) - curves can be distinguished by noting that the curves corresponding to the 
fastest evolution have the largest initial e. Increasing the eccentricity can be seen to reduce the inspiral time by up to several orders of 
magnitude over the circular case. In contrast, increasing the inclination in Fig.|4]only reduces the inspiral time by a factor ~ O(l) over 
a coplanar orbit. Also note that r e < Tj for all integrations. 



than the circularisation timescale (t, > r e ) unless a > 6. 
This was based on exponential decay estimates for small e 
and i, but nevertheless holds for the systems integrated in 
this work, as can be seen from Fig. [5] This makes intuitive 
sense, since a ~ 0(1) means that spin angular momentum 
is important, and circularisation involves only a property of 
the orbit, whereas alignment involves both the spin and the 
orbit. For typical HJs a < 6; we therefore expect r e < n, 
especially when the eccentricity damping effect of the tide in 
the planet is taken into account, which further enhances this 
inequality. The tide in the planet is completely negligible in 
changing the stellar obliquity since Ii <C h ~ 7712 a 2 . 



This means that if an orbit is initially inclined and ec- 
centric as a result of planet-planet scattering or Kozai mi- 
gration into a short-period orbit, we would expect the orbit 
to become circular before it aligns with the spin of the star. 
We should therefore observe fewer eccentric orbits than in- 
clined orbits, if those systems start with a umform distri- 
bution in e,i-space. This also means that if we observe a 
close-in planet on a circular orbit with non-zero i, we can- 
not rule out a non- negligible eccentricity in the past. 



8 APPLICATION TO OBSERVED SYSTEMS 

8.1 An explanation of the misaligned spin and 
orbit of XO-3 b 

The only system cu rrently observed with a spin-orbit mis- 
alignmenl^j is XO-3 ( |Hebrard et al. 



2008 



Johns-Krull et al. 



|2008|), which has a sky-projected spin-orbit misalignment 
angle of A ~ 70° ± 15°. This system has a very massive 
77i2 = 12.5 Mj planet on a moderately eccentric e = 0.29, 
P — 3.2 d orbit around an F-type star of mass 7711 = 1.3 Mq. 
Its age is estimated to be r* ~ (2.4 — 3.1) Gyr. Note that 
even if the star is rotating near breakup velocity (P* ~ 1 d), 
the planet is still subject to tidal inspiral, since P* > P cos i 
(where we henceforth assume i = X, which may slightly un- 
derestimate i). If we assume that the angle of inclination of 
the stellar equator to the plane of the sky is ~ 90°, then 
P* = 3.3 d ~ P i.e. fi ~ n. 



3 This spin-orbit misalignment has been confirmed since the sub- 
mission of the present paper by |Winn et al.||2009[ |, who find that 
the sky-projected spin-orbit misalignment angle A = 37.3°±3.7°, 
which is significantly smaller than that found by Hcbrard et al. 
|2008i. Nevertheless, our conclusions below should still be valid. 
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Figure 6. Tidal evolution of XO-3 b taking current values for 
the orbital properties of the system, except that cos i = 90° (not 
unreasonable since this roughly corresponds to the upper limit on 
A, which in any case gives a lower bound on i). Magnetic braking 
is included with 7 = 0.1, and Qi/n = 2 initially (results do not 
depend strongly on this choice). From the top and middle plots 
we require Q[ ^ 10 10 for the planet to survive for several Gyr, 
and maintain its high inclination. From the bottom plot we see 
that if Q' x ^ 10 10 , we require Q' 2 10 8 to maintain the current 
eccentricity for a few Gyr. Tidal dissipation in both the planet and 
star must therefore be weak to explain the current configuration 
of the system. 



Hebrard et al. ( 2008 1 quote a stellar spin-orbit align- 
ment timescale of ~ 10 yr for this system, but we find 
that this is in error by ~ 10 5 . We believe that the reason 
for this discrepancy is that their estimate was based on as- 
suming that the spin-orbit alignment time for XO-3 b is the 
same as for HD17156 b (Narita et al 



2008 Cochran et al. 



2008 1, which is a less massive planet on a much wider orbit. 
We find r, ~ 30 Myr (using the expression in assuming 
Qi — 10 6 to align the whole star with the orbit. Circularisa- 
tion time of r e ~ 10 Myr and the inspiral time is estimated 
to be r a ~ 16 Myr from simple estimates. 

Integrations for this system are given in Fig. [6] for a va- 
riety of stellar and planetary Q 1 values. These integrations 
again highlight the importance of considering coupled evolu- 
tion of the orbital and rotational elements, since timescales 
for tidal evolution can be quite different from the simple es- 
timates given above. Indeed, the actual spin-orbit alignment 
time from integrating the full set of equations is about an 
order of magnitude smaller than that from the simple decay 
estimate, due to the semi-major axis evolution. 

For the cases considered, the system can only survive 
and remain with its current inclination for ~ 3 Gyr if 
Q'x ^ 10 10 . An explanation for the survival and remnant 
orbital inclination of XO-3 b could therefore be the ineffi- 
ciency tidal dissipation in the host star. The host star is 
an F-star of mass mi = 1.3 ± 0.2Mq, so it will contain a 
small convective core and a very thin outer convection zone 
(OCZ) separated by a radiative zone. Dissipation in the con- 
vective core will only weakly affect the tide, and dissipation 
in the radiation zone will also be weak (Zahn 2008). This is 
because internal inertia-gravity waves excited at the inter- 
face between convective and radiative regions cannot reach 
the photosphere, where they can damp efficiently, as sup- 
posed for high-mass stars. In addition, nonlinear effects due 
to geometrical concentration of the waves in the centre of 
the star cannot occur because the waves will reflect from 
the outer boundary of the inner convection zone well before 
they become nonlinear ( jGoodman fc Dickson|1998| OL07). 

We expect that most dissipation occurs in the OCZ of 
the star. A calculation of Q' for the dissipation of the equilib- 
rium tide and inertial modes, in the thin OCZ, using a stellar 
modeQ appropriate for this star (note that the metallicity 
of the star is subsolar, with Z ~ 0.01), was performed (see 
OL07; OL04 for details of the numerical method). It must be 
noted that these calculations involve uncertainties regarding 
the effective viscosity of turbulent convection, though the 
general trends in the results below (and in the next section) 
are likely to be quite robust. The results for the m = 0, 1, 2 
components of the tide are plotted in Fig. [7] In Appendix 
[B| we show that a combination of the m — 0, 1,2 compo- 
nents of the tide are relevant for spin-orbit alignment and 
inspiral, and so we must calculate Q' for all components of 
the I = 2 (quadrupolar) tide. The relevant tidal frequencies, 
assuming Q ~ n currently, would be those of integer ui/fl. 
However, since the angle of inclination of the stellar equato- 
rial plane to the plane of the sky has not been determined, 
the relevant tidal frequencies cannot be calculated with any 
certainty. Nevertheless, Q' ^ 10 10 for most tidal frequencies 
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Figure 7. Tidal Q'-factor as a function of the ratio of tidal fre- 
quency to spin frequency lu/Q, from dissipation of the 1 = 2, 
m = 0, 1, 2 components respectively, of the equilibrium tide and 
dissipation of inertial modes in the OCZ of an F-type star (see 
OL07 for details of this calculation). This used a stellar model ap- 
propriate for XO-3 to model the convection zone. The dashed lines 
represent the effect of omitting the Coriolis force, and therefore 
inhibiting inertial waves. The prominent features in each figure 
(which occur for ui/Cl = —1 and —1/6 for m = 1 and ui/Q = — 1/3 
for m = 2) are Rossby waves, which are probably excited for tidal 
frequencies not relevant for the XO-3 system. For most tidal fre- 
quencies Q' 10 10 , which could explain the survival and remnant 

r^rtiifnl inrliTi n tinn nf ^fO_3 Vi 



for the host star XO-3. This can explain the survival and 
remnant inclination of XO-3 b, since both r a and t; are now 
much longer than the age of the system. 

In addition, the remnant eccentricity could be main- 
tained due to weak damping of the tide in the star for the 
same reasons. However, we must also explain the inefficient 
damping of the tide in the planet, if indeed the reason for 
the eccentricity is that r e > t*. The planet in this system is 
massive, and may be a low-mass brown dwarf. If it formed 
without a core, then the dissipation of inertial modes may 
be reduced if they are able to form global modes, as found 
in OL04. 

An alternative explanation for the survival of XO-3 b 
involving a larger initial semi-major axis for the planet is 
also possible, but this would require significant tidal migra- 
tion to bring the planet to its current location. This would 
require a much lower stellar Q' than our calculations predict, 
which as discussed in §3. 1| would imply a very short inspiral 
time for the planets on the tightest orbits, if their host stars 
have similar Q' values. We discuss this further with regards 
to the host stars WASP- 12 and OGLE-TR-56 in 30 below. 



8.2 Tidal dissipation in F— stars 

Our calculation of Q' for XO-3, indicates that it would be 
worthwhile to study the range of Q' expected for F-stars, 
and how these may differ from solar-type stars. We per- 
formed calculations of tidal dissipation in the OCZs of a 
variety of F-stars between the masses of 1.2 — 1.6 Mq, us- 
ing the numerical method of OL07, and the stellar models 
of EZ-Evolution. These stars contain convective cores sur- 
rounded by a radiative zone and an OCZ, and for the reasons 
mentioned in the previous section, we expect that tidal dis- 
sipation in the OCZ will dominate the dissipation. We con- 
sider the range of OCZ properties in F-stars, as the mass 
and metallicity are varied. Below we present a selection of il- 
lustrative examples, which represent the range of properties 
expected for F-stars. 

Fig. [8] shows Q' as a function of tidal frequency for two 
F-stars, with different OCZ properties. Fig. |8 (a)"] shows that 
tidal dissipation in more massive F-stars, with very thin 
OCZs, is found to be extremely weak, with Q' ^ 10 12 for 
most tidal frequencies. This implies that tidal dissipation in 
such stars is probably negligible in contributing to the spin- 



orbit evolution of HJs. Fig. 8(b) shows that tidal dissipation 



in lower mass F-stars with thicker convective envelopes, is 
similar to but slightly weaker than, that for solar-type stars 
(OL07), since Q' ~ 10 8 for most tidal frequencies in the 
range \Cj\ < 2|fl|. 

A model with 1.3 Mq star similar to XO-3, except that 
we choose supersolar metallicity (Z = 0.03), is found to 
have a similar OCZ to the 1.2 Mq model in Fig. [8(b)| and 
has a very similar Q' . The metallicity of the star affects the 
thickness of the OCZ, and therefore the efficiency of tidal 
disspation. 

From these examples, it is clear that assuming a single 
Q' applies for all stars is probably incorrect. Even within 
the mass range of F-stars there is considerable variation in 
Q' ~ 10 s - 10 12 in our calculations, primarily as a result 
of the variation in the mass fraction contained in the OCZ. 
There are also differences between G and F stars, due to dif- 
ferences in internal structure - most notably the presence of 
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Figure 8. Tidal Q'-factor as a function of the ratio of tidal fre- 
quency to spin frequency ui/Q, from dissipation of the I = m = 2 
component (we find the behaviour of the m = and 1 components 
is similar in magnitude) of the equilibrium tide, and dissipation of 
inertial modes, in the OCZ of two F-type star models. The dashed 
lines represent the effect of omitting the Coriolis force, and there- 
fore inhibiting inertial waves. The top figure shows a high— mass 
1.5 Mq F-star, with solar metallicity and an age of 0.7 Gyr, as- 
sumed to be spinning with a period of 3 d. This has a very thin 
OCZ, and tidal dissipation is extremely weak. The bottom figure 
shows a low— mass 1.2 Mq F-star, with solar metallicity and an 
age of 1.0 Gyr, assumed to be spinning with a period of 3 d. Tidal 
dissipation is much stronger in this star, though still smaller than 
that in Fig. 6 from OL07 for a 1.0 Mq G-star. Together with 
Fig.[7]this shows that Q' can vary considerably between different 
stars, even within the range of F-stars. 



a radiative core in a G-type star may lead to enhanced dis- 
sipation by nonlinear effects. Lower mass stars, and those 
with higher metallicity, tend to have thicker OCZs than 
higher mass, low-metallicity stars. In addition, higher mass 
stars are more centrally condensed, so the mass fraction in 
their outer regions will be lower. This results in low-mass, 
high-met allicity stars having lower Q' then high-mass, low- 
metallicity stars. In addition, OL07 found that the spin pe- 
riod of the star also affects Q' . 



8.3 The survival of WASP-12 b and OGLE-TR-56 
b 

The results of the last two sections also allow us to pro- 
pose an explanation for the survival of the planets on the 
tightest orbits, s uch as WASP-1 2 b ( |Hebb et aL]|2008| ) and 
OGLE-TR-56 b ( jSasselovpoM] ). Taking the current values 
for the stellar properties of the host stars in both systems, 
we find that they are both likely to have a similar internal 
structure to those discussed above, with a convective core 
present, which will prevent internal inertia-gravity waves 
from reaching the centre of each star. In addition, these stars 
are slowly rotating, and so the relevant tidal frequencies are 
likely to be outside the range of inertial waves > 2|n|), 
which would imply that such waves are not excited by tidal 
forcing. This means that Q' is likely to be much larger than 
that in the models discussed above (in Fig. [8|, particularly 
since reducing the spin frequency is found to increase Q' 
for a given ratio u/Q (OL07). Therefore it is likely that the 
relevant Q 1 > 10 10 , which would imply that tidally induced 
inspiral will not occur within the age of the system. This 
means that weak dissipation in the star could potentially 
explain the survival of both of these planets. 



9 CONCLUSIONS 

In this paper we have investigated the long-term tidal evolu- 
tion of close-in extrasolar planets. We studied the effects of 
magnetic braking on tidal evolution in a simplified system, 
and then performed numerical integrations for a variety of 
HJ systems, with particular emphasis on inclination evolu- 
tion. We now summarise the main results of this work. 

Magnetic braking moves the corotation radius of the 
star outwards such that any close-in planets will eventually 
orbit inside corotation (ficosi < n), and be subject to or- 
bital decay due to tides. It makes no sense to refer to any 
close-in planet as being tidally evolved, since magnetic brak- 
ing removes the only equilibrium state accessible through 
tidal friction during the expected stellar lifetime. 

Magnetic braking is found to be important for the tidal 
evolution of HJs unless Qcosi <C n, in which case the or- 
bit is already well inside the orbit-projected corotation ra- 
dius. Highly inclined (especially retrograde) orbits tend to 
be affected less by magnetic braking, since this condition is 
often satisfied regardless of the spin rate. Magnetic braking 
is particularly important when it comes to interpreting the 
tidal evolution of observed systems from formation to the 
present day, since the star may have been rotating much 
more rapidly in the past. Nevertheless, it is probably not a 
bad approximation to neglect magnetic braking for calculat- 
ing the future tidal evolution of most observed HJs, if the 
star has already spun down so that Q <C n. 

Combining our results with those of J ackson et al.| 



12008a), we find that coupled evolution of the orbital and 



rotational elements is essential to accurately model tidal evo- 
lution, and can result in much faster evolution than simple 
timescale estimates predict. This is especially true for highly 
eccentric orbits, for which the associated semi-major axis 
evolution increases the rate of stellar spin-orbit alignment 
by up to several orders of magnitude. 

We find that the true timescale for stellar spin-orbit 
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alignment is comparable to the inspiral time for HJs, there- 
fore the orbits of most close-in planets have probably not 
aligned with the spin of the star. Observed inclinations are 
likely to be a relic of the migration process. This means that 
RM observations of transiting planets can potentially distin- 
guish between migration caused by planet-planet scattering 
or Kozai oscillations combined with tidal dissipation in the 
star, and that produced by tidal interaction with the gas 
disc. If the majority of candidates are found with A ~ 0, 
this strongly disfavours planet-planet scattering or Kozai 
migration, since they are expected to produce significantly 
inclined orbits, and we have found that tides are unlikely to 
have aligned orbits without causing inspiral. Alternatively, 
if systems are found with significantly nonzero A, then some 
planet-planet scattering or Kozai migration could have oc- 
curred to produce these orbital inclinations. We strongly en- 
courage future observations of the RM effect for transiting 
planets. 

For most HJs, tides tends to circularise the planet's 
orbit before spin-orbit alignment or inspiral occurs. If an 
orbit is initially inclined and eccentric, as a result of planet- 
planet scattering or Kozai migration into a short-period or- 
bit, then we would expect the orbit to become circular before 
it aligns. Therefore, if we observe a planet on an inclined, 
circular orbit, we cannot rule out a non-negligible eccentric- 
ity in the past. This means that we should observe fewer 
eccentric orbits than inclined orbits if those systems start 
with a uniform distribution of points in (e,i)-space, due to 
tidal friction. This should be considered when comparing the 
observed (e, i) distribution with those predicted from theo- 
retical work on Kozai migration or planet-planet scattering, 



before we can further constrain these theories (Fabrycky & 
Tremaine|2007l |Juric fc Tremaine|2008 1 



The misaligned spin and orbit of the XO-3 system could 
potentially be explained in terms of inefficient tidal dissipa- 
tion inside the host star. The required stellar Qi > 10 10 
required for the survival and remnant misalignment is pre- 
dicted from theoretical calculations of tidal dissipation in 
the OCZ of an F-star. In addition, the remnant eccentricity 
poses constraints on the planetary Q' 2 > 10 8 , in the absence 
of perturbing forces that could excite the eccentricity. 

The stellar Q' has been shown to vary widely between 
different stars, and even within the range of F-stars, from 
10 8 - 10 12 . There are also differences between Q' for the F- 
star models discussed here, and the solar-type star models 
discussed in OL07. This implies that assuming a single value 
of Q' applies to all exoplanet host stars is probably incor- 
rect. Q' has been found to vary with stellar mass, metallic- 
ity, spin period, as well as tidal frequency. The presence of 
a convective core is likely to be important, in that it acts 
to prevent internal inertia-gravity waves, which are excited 
at the interface between the convective and radiative zones, 
from reaching the centre of the star, where nonlinear effects 
could enhance their dissipation. This may explain the sur- 
vival of some of the planets on the tightest orbits, such as 
WASP-12 b and OGLE-TR-56 b, whose host stars are likely 
to have convective cores. 

It is clear that much more work is required to study the 
mechanisms of tidal dissipation in rotating stars. It would 
also be useful to perform a detailed study into the possibility 
of core-envelope decoupling in such stars, since this may 
have implications for the survival and remnant inclinations 



of close-in planets. Future observations of transiting planets 
will constrain these theories, hopefully leading to a better 
understanding of the mechanisms at work in the formation 
and evolution of planetary systems. 
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APPENDIX A: TIDAL EVOLUTION EQUATIONS 



Here we present the orbit-averaged tidal evolution equations (Eggleton et al. 1998 Egglcton & Kiseleva-Eggleton 2001 
Mardling & Lin 2002 1 used in this work. These are derived starting with the equation of relative motion of a planet of mass 
m,2 and its host star of mass mi 



d^r 

df 2 



Gmis 



r + f, 



(Al) 



where r is the separation, mi2 = mi + m2 , and f represents a perturbing acceleration. In Egglcton & Kisclcva-Egglcton ( 2001 ) , 
f contains many contributions: the effects of general relativity, perturbing accelerations due to other planets and tidal and 
spin distortions of the star and planet, as well as tidal friction. Here we solely consider the perturbing effects of tidal friction 
and set f = f t f = f t f + f t 2 f , where 



-3ri,2fci,2W 



[ m2 > 1 
\mi } 2 



Rl,2 

a 



[3(r ■ r)f + (f x r — r fl^a) X f] . 



(A2) 



Here n = y \ii 12 i s the orbital mean motion, ki t 2 are the second-order potential Love numbers for the star and planet 
respectively (which is twice the apsidal motion constant), and ri,2 is the effective tidal lag time for each body. This form of 
the dissipative force of tidal friction is that derived under the assumption of a constant lag time in the equilibrium tide model 
(Egglcton et al.|1998 1 . In this model, we assume that the body quasi-hydrostatically adjusts to the perturbing potential of its 
companion, but delayed by some small lag time (n for the star, T2 for the planet) that is proportional to the dissipation. Thus 
for each body, Q = ^ is assumed to be inversely proportional to the tidal frequency u), so that the lag time r is independent 
of tidal frequency, and is therefore the same for all components of the tide. This model is beneficial because the resulting 
secular evolution equations can treat arbitrary orbital eccentricities and stellar and planetary obliquities. 

In the resulting equations we have chosen to parametrise the efficiency of tidal dissipation in each body by redefining 
Q = — , and adopt a constant Q which does not change during the evolution i.e. we assume effectively that the lag time scales 



with the orbital period (then Eq. A2 matches Eq. 4 in Mar dling & Lin 2002 where this assumption was not made explicit). 
We also introduce the definition Q = as discussed in £3.1 This allows us to discuss a Q' values" for particular bodies, 
which do not change as the orbital and rotational elements vary. For the purposes of this paper we define Q' = 2 , 3 . This 
is equivalent to assuming that the relevant tidal frequency ill = n. Note that this may not give identical numerical factors 
in the resulting equations to other formulations of tidal friction (e.g. Goldreich & Soter 1966; Zahn |l977 Hut|1981~ l, but we 
feel that this is the best way to study the general effects of tidal friction, given our uncertainties in the value of Q , and its 
dependence on Cj, for realistic giant planets and stars. 

The secular evolution of the orbital elements is calculated via the specific orbital angular momentum vector h = rxr = 
na 2 \/l — e 2 h, and the eccentricity vector e. The eccentricity vector has the magnitude of the eccentricity, and points in the 
direction of periastron, and is given by 



r x h 

Gmi2 



(A3) 



The planet's specific orbital angular momentum h changes at a rate 

dh 

~dt 



r x (f^ + f t 2 f ) 



(A4) 



with a corresponding rate of angular momentum transfer between the orbit and spin of each body given by 



Jl,2 = /l, 2^1,2 = x f tf > 



(A5) 



since total angular momentum is conserved, and where n 
evolves as 



is the reduced mass of the system. The eccentricity vector 



rfe = [2(f tf ■ r)r - (r ■ r)f tf - (f tf ■ r)r] 
dt Gm\2 



(AQ) 



Eqs. |A4| and |A6| are time-averaged over the orbit of the planet, and the resulting differential equations are given below. 
Numerical integration of these equations gives the secular evolution of the orbital elements. Note that we have written these 
equations so that they are regular at e = 0, unlike those in Eggleton & Kiseleva-Eggleton (20011 and Mardling & Lin (20021 



- we eliminate reference to the basis vectors chosen in their representation, since e is undefined for a circular orbit, whereas 
e is well defined and equal to zero. 
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(A9) 
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We also need to define an inverse tidal friction timescale t^ 1 for each body (here for body 1, change 1 — > 2 to get the 
corresponding expression for body 2), and the functions of the eccentricity (first derived in a similar form by Hut|1981 |. 
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(A21) 

(A22) 
(A23) 

In the absence of magnetic braking (o; m b) the total angular momentum is conserved (4r = 0). With the inclusion of magnetic 
braking, the total angular momentum of the system decreases at a rate 

= — Quh + Jifii + J 2 fl 2 ) = Jitjmb = -amblxfifOi. (A24) 

Tidal dissipation and magnetic braking result in a loss of energy from the system. The total rate of energy dissipated can 
be expressed as follows: 

1 Gm\m,2 d 
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where the final term in the last expression is the energy dissipated by magnetic braking and 
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are the tidal heating rates in the star and planet, respectively. It has been proposed that H2 can result in planetary inflation, 



On the tidal evolution of Hot Jupiters on inclined orbits 17 



if the energy is deposited deep in the interior of the planet ( Bodenheimer et al.||200i Jackson et al.|2008b \ . Note that these 
expressions approach zero as e — > 0, fli,2 ■ e — > and f2i,2 • h — > n i.e. tidal dissipation continues until the orbit becomes 
circular, and the spin of each body becomes coplanar and synchronous with the orbit. 



APPENDIX B: TIDAL POTENTIAL EXPANDED TO ARBITRARY ORDER IN STELLAR 
OBLIQUITY 

In this section we expand the tidal potential into its separate Fourier components, following the approach of OL04. We aim 
to study which tidal frequencies are relevant for highly inclined orbits, such as that of XO-3 b. We consider two bodies in 
mutual Keplerian orbit with semi-major axis a and mean motion n. Adopt a coordinate system with origin at the centre of 
body 1, and let this represent the star, and body 2 the planet. The tidal potential experienced at an arbitrary point P in 
body 1 is the nontrivial term of lowest order in r, 

* = ^[«V-3(R.r)*], (Bl) 

where position vector of the point P in body 1 is r, and the position vector of the centre of mass of body 2 is R(i). Body 2 is 
treated as a point mass, of mass mi. We consider an inclined, circular orbit. Without loss of generality, we consider body 2 
to orbit in a plane inclined to the (x, j/)-plane by an angle i, so that its Cartesian coordinates are 



R = a (cos nt cos i, sin nt, cos nt sin i) , 

while the Cartesian coordinates of a point P are 

r = r (sin 9 cos <f>, sin 9 sin <f>, cos 9) , 

where (r, 9, <j>) are the usual spherical polar coordinates. Then 
* = ^ m , 2 r 2 [l — 3 (cos nt cos i sin 9 cos (f> + sinnt sin 9 sin <j> + cos nt sin i cos 9) 2 ] 
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where ^ m ^ I and I g Z + , denote an associated Legendre polynomial normalised such that 
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The tidal potential correct to arbitrary order in the stellar obliquity i can be expanded as a series of rigidly rotating solid 
spherical harmonics of second degree, 
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where the azimuthal order n%j, frequency uij, tidal (Doppler-shifted) forcing frequency (In- 
dependent amplitude Aj(i) of each component are 
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cos = — 2n, a>8 = — 2n — 2fi, As = — 
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(B15) 



There are eight components of the tide that contribute for arbitrary stellar obliquity. For small inclination, the terms that are 
of first order in the obliquity are for m = 1 and are the j = 4 and j = 5 components above (not the component having tidal 
frequency u> = n — Q, as incorrectly stated in OL04). For a coplanar orbit, only j = 1 and j = 2 components above contribute, 
and these reduce to the first two components in OL04. 

For an orbit at i ~ 90°, all components except j = 2 and j = 4 contribute to the tidal force. Thus, for XO-3, it is important 
to calculate the resulting Q' for the m = 0, m = 1 and m = 2 components of the tide. The relevant tidal frequencies for this 
system cannot be calculated with any certainty, since the stellar spin period has not been accurately determined. If we assume 
that the angle of inclination of the stellar equatorial plane to the plane of the sky is ~ 90°, then we have O ~ n currently. 
The relevant tidal frequencies would then be to = 0, ±fi, ±2fi, ±3fi, ±4f2. 
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